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(57) Abstract 

A novel method for estimating the 
end-systolic pressure-volume relationship 
(ESPVR) is presented. The method provides 
for accurate estimation of the ESPVR from 
a single beat of the cardiac cycle. The 
method is based on normalized human 
time varying elastance curves [En^n)]. The 
ESPVR is estimated from one beat using 
PV data measured at normalized time tN and 
end-systole (t m ax) to predict intercept: Vn(SB) 
« tPNONJxVamax) - V(t N )xE N (t N )]/tPN(tN) - 
EnOn)] and slope: E max (SB) = P(traax)/[V(tma*) 

- Vo(SB)]. The present invention provides 
a method for ESPVR estimation which is 
non-invasive and particularly applicable for 
bedside applications. 
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ASSESSING CARDIAC CONTRACTILITY 
AND CARDIOVASCULAR INTERACTION 

BACKGROUND OF THE INVENTION 

5 

1 . Field of the Invention 

The present invention relates generally to the field of cardiovascular 
assessment and performance. More specifically, the present invention 
relates to the field of non-invasive assessment of cardiac function based on 
10 pressure-volume relation analysis. 

2. Description of the Related Art 

The pressure-volume framework for assessing cardiac function, 

15 specifically the use of an end-systolic pressure-volume relationship (ESPVR) 
to assess cardiac contractile function is well-known. This analytic method is 
widely regarded as among the best means for defining cardiac chamber 
performance, and its interaction with the vascular (venous and arterial) 
systems. It can be used to both predict the behavior of the cardiovascular 

20 system during physiologic changes (such as sudden increases or decreases 
in blood pressure or circulating blood volume, and the influence of 
pharmacologic agents typically used to treat patients with heart disease. The 
major impediment to practical implementation of this analysis for clinical 
medicine is the difficulty in making all of the necessary measurements, and 

25 the requirement for invasive data. Recent advances in pressure-monitoring 
and Doppler imaging technology have provided methods to estimate arterial 
vascular properties non-invasively, but ESPVR characterization has remained 
much more difficult. A key contributor to this difficulty is the requirement of 
collecting data from many cardiac cycles for which the "loading" conditions 

30 (filling volumes, ejecting pressures) of the heart have been altered. 

There have been attempts to estimate the ESPVR from single-cardiac 
cycles, however, these too do not provide a reliable and easily implemented 
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method. The development of single-beat ESPVR estimation methods has 
been driven by the desire to simplify the loading procedures required to 
directly measure the ESPVR, and to avoid a need for continuous LV (left 
ventricular) volume data. In humans, the ESPVR is for all practical purposes 
5 a linear relation, which can be defined by a slope (termed either E max , or the 
end-systolic elastance, E es ), and a volume axis intercept (V 0 ). The slope is a 
measure of chamber stiffness, and can be used to index contractility changes, 
while the volume axis intercept is a measure of chamber geometry, and can 
also be used to index contractile performance changes. Two basic 

10 approaches have been employed to date. The simplest and most frequently 
used has been to ignore V 0 , and report the ratio of Pes/Ves (end-systolic 
pressure/volume ratio). This has shown to be inadequate since there can be 
marked variability in V 0 , with often large non-zero intercepts in patients with 
infarction or dilated cardiomyopathy. An alternative method uses abrupt 

15 occlusion of the aorta during isovolumic contraction to determine peak 

isovolumic pressure, LV volume during ejection is estimated by integrating 
aortic flow and synchronizing the result with LV pressure. This method yields 
two pressure-volume points, one for the isovolumic contraction and one for 
the resting ejecting contraction, and by linking the two points by a line, one 

20 can estimate the slope E max . A modification to this method was first 

suggested by Sunagawa et al, ". . .Instead of actually occluding the Aorta to 
determine isovolumic pressure and volume, an approach that could never be 
used clinically, these investigators provided a method to estimate this 
pressure-volume (PV) data by mathematical curve fitting of LV data measured 

25 during isovolumic contraction and relaxation phases of a steady-state ejecting 
beat." The fit predicts the peak pressure that would be generated in the 
absence of ejection, and by drawing a line between this theoretically derived 
"isovolumic" pressure-volume point and the measured end-systolic PV point 
from the same ejecting beat it yields the ESPVR estimate. This method has 

30 been tested in normal animals and in humans. While providing adequate 

results, the method required measurements of high fidelity LV pressure and is 
based on an arbitrary curve fit. Current non-invasive methods for assessing 
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heart contractile function and the interaction between the heart and the 
arterial system are not easily implemented in a clinical environment. 

SUMMARY OF THE INVENTION 

It is an object of the present invention to provide a method for 
accurately predicting ESPVR from LV PV data that is easily extracted from a 
single heart beat. 

It is a further object of the present invention to provide a method for 
accurately determining ESPVR which may be implemented non-invasively at 
the bedside. 

It is another object of the present invention to provide a method to link 
the ESPVR predictions with non-invasive imaging and pressure monitoring 
technologies to generate clinical displays that provide a treating physician 
with information about current cardiovascular status and the likely effects of 
pharmacologic therapies on this status in patients with heart disease. 

The end-systolic pressure-volume relation (ESPVR) is a valuable 
measure of ventricular systolic function for both clinical and experimental 
evaluations providing a very useful measure of contractile function. However, 
until now, in order to utilize ESPVR, it has been necessary to acquire data 
from multiple cardiac cycles at varying loads. This has limited our ability to 
utilize ESPVR. Acute ESPVR shifts primarily indicate inotropic change. The 
slope (E max ) measures end-systolic chamber stiffness, an important 
determinant of the influence of vascular loading on systolic pressure and flow 
generation. The ESPVR evolved from a linear time-varying elastance [E(t)] 
model of contraction, in which the left ventricle (LV) was considered to be an 
elastic structure that stiffens with a predictable time-course during the cardiac 
cycle. 

The present invention utilizes a normalized time varying elastance 
model, E(t) curve (normalized both by amplitude and time to peak amplitude) 
to approximate the ESPVR from pressure and volume data relating to a single 
heart beat. Based upon collected data from patients reflecting a wide range 
of cardiovascular diseases as well as normal hearts, a normalized time- 
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varying elastance curve was derived. The normalized E(t) function was 
defined based upon the elastance curve. The time interval between the R- 
wave of the electrocardiogram (ventricular activation) and the maximal PA/ 
ratio (t^) was determined, and real-time coordinates transformed to 
5 normalize time by this interval. The ESPVR estimate is based on two 
assumptions: 1) that the E(t) model is linear, and 2) that a single and 
constant volume axis intercept applies for a given cardiac cycle. 

Based upon equations for a time-varying elastance, and specifically 
the ratio of instantaneous chamber elastance to maximal elastance measured 
10 at the end-systole, the ESPVR volume axis intercept was calculated. From 
this determination, the slope is calculated, providing the necessary numbers 
to estimate the ESPVR. 

BRIEF DESCRIPTION OF THE DRAWINGS 

15 Figure 1 is a schematic diagram explaining estimation of ESPVR slope 

and intercept. 

Figure 2 is a diagram of normalized elastance [E N (t N )] curves for 
various patient groups and/or test condition for the underlying data supporting 
the present invention. 
20 Figure 3 is a diagram illustrating group averaged normalized elastance 

curve [E N (t N )] from the underlying data supporting the present invention. 

Figure 4 shows comparisons between measured values based on 
underlying patient studies and values determined by the present invention. 

Figure 5 compares a single-beat ESPVR slope estimation using 
25 invasive arterial pressure data and ventricular pressure data. 

Figure 6A illustrates a first embodiment of the invention which 
measures pressure and volume using invasive methods. 

Figure 6B explains the flow of data and the data processing that is 
carried out according to the first embodiment of the invention. 
30 Figure 7 illustrates a second embodiment of the invention in which 

pressure and volume are measured using non-invasive methods. 
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Figure 7B explains the flow of data and the data processing that is 
carried out according to the second embodiment of the invention. 

Figure 8a is a diagram illustrating an application of the single-beat 
ESPVR estimation method of the present invention. 
5 Figure 8b is a diagram illustrating an application of the single-beat 

ESPVR estimation method of the present invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

It has been found that E N (t N ) curves are surprisingly similar among 

10 _ hearts, particularly during early contraction. Based on this observation, the 
present invention provides an arrangement (method and apparatus) for 
determining a good approximation of ESPVR using pressure, volume and 
timing data easily obtained from a single cardiac cycle. Our arrangement 
identifies both resting and acute changes in the ESPVR, is minimally 

1 5 sensitive to load, and yields a better estimate of directly measured relations 
than previously proposed techniques. Further, the present invention is 
adaptable to non-invasive analysis. As such, our invention makes it possible 
to more easily utilize ESPVR analysis in clinical medicine. The method has 
value not only for assessing acute and chronic changes in contractile function 

20 from therapeutic interventions, but also for predicting a patient's 
hemodynamic response to therapy. 

The claimed arrangement for single-beat ESPVR estimation is based 
on the underlying premise of the time-varying elastance model of cardiac 
contraction that there is consistent E N (t N ), independent of heart rate, loading 

25 conditions, or contractility, that is scaled by both amplitude (E,^) and time 

(tmax) fo r a 9 ive n ventricle. This further implies that the elastance relations are 
linear, and that they intersect at a common volume intercept. 

The underlying data which was collected to confirm the basic 
hypothesis for the present invention established that the E N (t N ) curves were 

30 consistent among patients despite markedly varying underlying diseases and 
conditions, and that a single-beat estimation that assumed ESPVR linearity 
and constant V 0 yields very reasonable predictions. Establishing that E N (t N ) is 
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consistent among hearts despite widely varying chronic disease states is 
what sets the present invention apart. 

The present invention is based upon a study of nearly one hundred 
patients having a variety of heart conditions, including normal hearts and 
5 varying forms of chronic heart disease. Using these patients a reference or 
"gold standard" ESPVR was derived from multiple variably preloaded cardiac 
cycles measured during transient obstruction of IVC (inferior vena cava) 
inflow. Referring now to Figure 1A, following the transient obstruction of the 
IVC, the blood return to the left-ventricle is decreased and a series of cardiac 

10 cycles with decreasing volume and pressure results. The ESPVR line 
connects the upper-left corners of the PV loops, which are the points of 
maximal pressure/(volume-V 0 ), where V 0 is the ESPVR volume-axis intercept. 
Following the determination of the ESPVR from multiple beats (MB), its 
parameters , the slope E max(MB)l and the volume axis intercept V 0(M b). are 

15 directly calculated. Figure 1B presents how pressure-volume data from a 
single beat (SB) loop, which is either measured from one cardiac cycle or 
averaged from several cardiac cycles during steady state, are used by the 
current invention to estimate the ESPVR parameters. Based on two point on 
the PV loop, at times t N and t max , the slope E max(SB) and intercept V 0(SB) are 

20 estimated as detailed below. 

To derive the normalized time-varying elastance curve, PV data from 3 
to 5 sequential steady-state beats were signal averaged to yield a PV loop. 
Time-varying elastance [E(t)] was defined as the instantaneous ratio of 
P(t)/[V(t)-Vo(MB)I. using V 0(MB) from prior multiple-beat ESPVR analysis. The 

25 maximal value of E(t) [E max(SB) ] and the time to achieve E max(SB) referenced from 
the R-wave of the electrocardiogram (t^). were both determined. The 
normalized E(t) function was then defined as: 

E N (t N ) = E(t N )/E max(SB) Eq. 1 

where t N = t/t^ Eq. 2 

30 To combine E N (t N ) relations, each curve was re-sampled at 200 equispaced 
intervals using linear interpolation, and the results averaged for each patient 
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and/or hemodynamic group. A total mean E N (t N ) curve was obtained by 
averaging the results. 

The present invention provides a method for single-beat ESPVR 
estimation, characterized by a linear slope E^g, and a volume-axis intercept 
5 V 0(SB) . Figure 1 B illustrates schematically the estimation method of the 

present invention. For each steady-state cardiac cycle, t max (the time between 
the R-wave and the maximal PA/ ratio) is first determined, and the time 
coordinates are transformed to normalized time t^ utilizing Eq. 2. PV data is 
then measured at times t N [P(t N ) and V(t N )] and t^ [P(t max ), Vft^)], using 

10 linear interpolation. The ESPVR estimate is based on two assumptions: 1) 
that the E(t) model is linear; and 2) that a single volume axis intercept (V 0(SB) ) 
applies for a given cardiac cycle. Using the definition for time-varying 
elastance E(t), chamber elastance at time t N is: 

E(t N ) = P(t N )/[V(t N )-V 0(SB) ] Eq.3 

1 5 and elastance at time t^ is: 

E(tmax) = Emax(SB) = P(UWJ - V 0(SB) ] Eq. 4 

As established, the ratio of E(t N )/E max{SB) equals E N (t N ). Therefore, combining 
Eqs. 1, 3 and 4 the volume-axis intercept is obtained as: 

20 V 0(SB) = [P N (t N )xV(t m J - V(t N )xE N (t N )] 

[P N (U - E N (t N )] Eq.5 

where P N (t N ) = PfW/Pft^). Once V 0(SB) is calculated, E max{SB) is determined 
from Eq. 4. 

25 The results of the underlying study confirm the method of the present 

invention. With respect to the normalized elastance curve, E N (t N ), Figure 2 
displays averaged E N (t N ) curves for patients with normal hearts (A), and those 
with hypertrophy, (B), dilated cardiomyopathy, (C), coronary disease (D), and 
post-infarction aneurysm (E). The solid line is the mean and the dashed lines 

30 show ± 1 Standard Deviation (SD). There is relatively little variance about the 
mean for each E N (t N ) curve in the various patient groups, particularly during 
the initial contraction phase (prior to the first inflection point). Furthermore, 
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the mean E N (t N ) curves for the different groups were very similar, despite 
marked differences in preload, afterload, contractility and chronic disease. 
This is shown in Figure 2H in which the curves are superimposed. E N (t N ) 
curves measured before or after positive inotropic stimulation (Figure 2F, E max 
5 from 1 .5 ± 0.5 to 2.7 ± 0.9) or artrial tachycardia pacing (Figure 2G, HR at 76 
±10, 99 ±9 and 130 ±12 

min" 1 ) were virtually superimposable. The phase with the least variance is 
early contraction (isovolumic phase), while that with the largest occurs during 
late relaxation. The net averaged E N (t N ) ± 1 standard deviation (SD) curve is 

10 shown in Figure 3. 

The results of the underlying study also establish the consistency of 
E N (t N ) supporting the method of the present invention for single-beat ESPVR 
estimation. While Eqs. 4 and 5 yield V 0(SB) and E max(SB) estimates for any time 
t N , it is preferable to select a t N when the physiologic variance in E N (t N ) is 

15 minimal. From Figure 3, this occurs during early contraction. In addition, Eq. 
5 reveals V 0(SB ) to be a hyperbolic function of E N (t N ) with a vertical asymptote 
at E N (t N ) = P N (t N ). If t N (and E N (t N )) is selected near this value, then small 
errors will result in greater instability of the V 0(SB) estimate. The optimal t N is 
one that minimizes both physiologic and mathematical variances. It is 

20 preferred that t N lie between 0.25 and 0.35. Figure 4 displays results 

comparing directly measured values for the slope and intercept of the ESPVR 
with those values predicted by the current invention. Figure 4A compares 
E max and V 0 estimated with the new single beat method (E max(SB) and V 0(SB) ) 
with the "gold standard" estimation based on multiple beats (E max(MB) and 

25 V 0(MB) ). The solid line is the linear regression line , and the dashed lines 

represent the 95% prediction intervals. Both sets of parameters were highly 
correlated, with regression line given by E^sa, = 1.01 X E max(MB) + 0.03 , 
correlation coefficient r=0.9, standard error of estimate is 0.42. Figure 4B 
presents the difference between the single beat and multiple beat estimations 

30 as function of the multiple beat estimations, and shows minimal dependence 
between the difference and the "gold reference" values of E^ and V 0 . The 
new method provides good predictive value, with a small standard error of the 
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estimate, and regression slope and intercept that are very close to 1.0 and 
0.0 respectively. 

The ability to provide accurate non-invasive ESPVR estimation is 
particularly desirable. Clinical evaluation of ventricular function and 
5 ventricular-arterial interaction is particularly useful when employed as part of 
the chronic assessment of heart disease. The ability to make measurements 
non-invasively is valuable, since repeated evaluations over time generally 
provide the most insight, and these measurements are impractical if invasive 
data is required. Since the method of the present invention provides ESPVR 
10 estimation which incorporates data near the end of isovolumic contraction, it 
is particularly suitable for non-invasive assessment. Echocardiography or 
nuclear ventriculography can provide end-systolic and end-diastolic chamber 
volumes, while arterial cuff pressures provide diastolic and systolic pressures. 
Timing of the onset of ejection and end ejection can be obtained non- 
15 invasively by Doppler cardiography, by Phonocardiography, or other methods. 

The ability to incorporate parameters which may be measured non- 
invasively provides a significant advantage to the present invention. E max(SB) 
may be estimated using Aortic pressure data. By selecting t N to occur at the 
onset of cardiac ejection then: V(t N ) = EDV, and P(t N ) = AJ^,*, the Aortic 
20 diastolic pressure. The end-systolic pressure, or P(t max ) could be estimated 
from the arterial systolic pressure by: 

P(t m J = Pes = (0.9xA 0 P SYS ) Eq.6 
Therefore, from arterial systolic and diastolic pressures, and values for ESV 
and EDV, each being measurable non-invasively, V^ can be estimated by: 
25 V 0(SB) ^ OP = [A o P DlA xESV/(0.9xA o P SYS )-EDVx E N (t N )] 

[A o P DlA /(0.9xA o P SYS )- E N (t N )] Eq.7 

and 

E ma x ( sB)-AOP = (0.9xA o P SYS )/(ESV-V 0(SB) ) Eq. 8 

For an analysis based on arterial pressures, t^ is calculated as the 
30 time between the R-wave of the electrocardiogram and the time of end- 
systole, which can be approximated by a time when arterial pressure declined 
by 10% of peak or by a time of end of ejection. Further, t N is the normalized 
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time at the upstroke of arterial pressure (i.e. diastolic pressure) and E N (t N ) is 
determined at this time based on the previously derived averaged E N (t N ) 
curve, and the value is utilized in Eq. 7. Figure 5 illustrates that with accurate 
measurements and timing, this substitution yields excellent agreement. 
5 The invention can be used for the estimation of the ESPVR from either 

invasive or non-invasive measurements. One preferred embodiment is for 
use in the catheterization laboratory, where invasive pressure-volume 
measurements are generally recorded. Referring now to Figure 6A, the 
apparatus includes a high fidelity invasive pressure measurement system (10) 
10 which provides a continuous measurement of the left-ventricular blood 

pressure but preferably provides also the aortic blood pressure. Catheters 
with dual micromanometer sensors, or single sensors that provide these 
signals are readily available (e.g. catheter model SPC-320, Millar Instruments, 
USA), as are the transducers and interface electronics to pass these signals 
15 to all current commercial catheterization laboratory hemodynamic recording 
systems (e.g. Q-Cath, Quinton Instrument Company, USA). Secondly, the 
apparatus involves an invasive volume measurement system (12) (this can be 
based on intraventricular conductance measurement system like Sigma V, 
CardioDynamics, the Netherlands, which gets the conductance 
20 measurements through conductance catheter, e.g. SSD-767. Millar 

Instruments, USA), that yields a continuous left ventricular volume signal. 
Third, is a standard electrocardiographic monitoring system (14) that provides 
a continuous analog ECG signal from a selected lead configuration, and is 
typically part of the commercial catheterization laboratory hemodynamic 
25 recording systems. 

The pressure (16), volume (18) and the electrocardiogram (20) analog 
signals are digitized in real-time by a standard A/D board (22) (e.g. AT-MIO- 
16E-10, National Instruments, USA), and are displayed in real-time as well as 
recorded to computer memory using a standard personal computer (24) (e.g. 
30 Pentium 166MHz PC, Gateway 2000, USA). The signals can be displayed 
using signal acquisition software such as provided by LabView 4 (National 
Instruments, USA) or by use of custom designed software. 
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Referring now to Figure 6B, the analysis of the digitized pressure (30), 
volume (32) and ECG (34) signals by the personal computer (24) is 
schematically presented and explained below. 

Pressure (30) and volume (32) signals are recorded during several 
5 cardiac cycles (typically 3-5 cycles), ideally under conditions of quiet end- 
expiration to minimize signal noise. The onset of each cardiac cycle is 
determined by identifying the R-wave (36) of the ECG (34) through standard 
triggering modes such as thresholding of the ECG signal or its derivative. 
Based on the identification of the R-wave (36), the digitized pressure (30) and 

10 volume (32) signals are sequenced into pressure_and volume signals of 

separate cardiac cycles. The pressure (38) and volume (40) signals from the 
separate cardiac cycles are group averaged to yield a mean pressure (42) 
and mean volume (44) representing the pressure and volume change over a 
cardiac cycle. From the average pressure (42) and average volume (44) 

1 5 data, the analysis program identifies the occurrence of end-ejection, based on 
well known physiologic features of these signals (for example, the end- 
ejection point can be identified from the volume signal as the point where 
there is no more volume decrease due to ejection, or it can be identified as 
the point where the aortic pressure and the ventricular pressure deviate 

20 following the closure of the aortic valve). The time interval between the 
occurrence of the R-wave (36) and end of ejection is denoted t^ (46). 

A default value or user-controlled value of the normalized time t N (48), 
which from experimental data should be within the range of 0.25 to 0.35, is 
used to calculate the absolute time T (50) which corresponds to the 

25 normalized time t N . The absolute time T is calculated by rearranging 

Equation 2 to the form: t = t N x t^ . The pressures (52) and volumes (54) 
corresponding to times T (50) and "t max " (46) are calculated by interpolation 
from the averaged pressure (42) and volume (44) data. The pressure at time 
T is normalized by the pressure at time t max and the normalized pressure is 

30 denoted P N (t N ) (60). The volume at absolute time T, which corresponds to 

normalized time ^ is denoted V(t N ). The volume at time t max is denoted V(t max ). 
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The normalized elastance function for time t N , namely E N (t N ) (62) is 
determined from the stored mean elastance function (64). For example, if the 
mean elastance function is stored in tabulated form, a linear interpolation can 
be used to find the E N (t N ) from the nearest tabulated points. If the mean 
5 elastance function is stored as a formula (e.g. polynomial or Fourier fit), the 
E N t( N ) can be calculated directly. 

Finally, the volume-axis intercept of the ESPVR is calculated by 
Equation 5, and the slope is calculated by Equation 4. The calculated slope 
and intercept of the ESPVR (66) are presented to the user as numerical 

10 values and by a graphic display as in Figure 4r 

A second preferred embodiment is for use with non-invasive 
measurements. Referring now to Figure 7A, the apparatus includes either a 
standard or a custom made non-invasive blood pressure monitor system (70) 
(e.g. Dynamap 1846SX, Critikon, USA ). This monitor provides systolic and 

15 diastolic arterial blood pressure. Other alternative pressure signals would be 
continuous pulse tracings obtained by carotid tonometry (SPC350, Millar 
Instruments, USA), or by transformed radial arterial pulse tracings obtained 
by radial tonometry (Jentow, Collin Medical Instrument Corporation, USA). 
The second component (72) is a standard echo-Doppler system (e.g. Sonos 

20 2500, Hewlett-Packard, USA). The commercially available echo-doppler 
systems provide clinically accepted estimates of ventricular end-diastolic 
volume (EDV), end-systolic volume (ESV) by standard analysis tools. An 
alternative option is to use a digital edge-detection algorithm, which is 
available as an option on some echo-doppler systems (e.g. Acoustic 

25 Quantification option, Sonos 2500, Hewlett-Packard, USA) to provide 
calibrated volume waveforms - in addition to the EDV and ESV values. 
Doppler flow measurement is used to derive stroke volume by time integral of 
the aortic velocity flow measurement multiplied by aortic cross sectional area, 
which is a standard clinical approach and is implemented as an analysis tool 

30 on most cardiologic echo-doppler systems. Lastly, the ejection fraction is 
derived from measurements using M-mode or directly from 2-dimensional 
imaging, both modes are standard on commercial echo-doppler systems. 
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From the ejection fraction (EF) and stroke volume (SV) , one calculates the 
end-diastolic (EDV) and end-systolic volumes (ESV) from: EDV = SV/EF; 
and ESV = SV(1-EF)/EF, based on the definitions of these parameters. 
These calculations are also provided as standard analysis tools by most 
5 cardiologic echo-doppler systems. Third, the system involves an 

electrocardiographic recorder or monitor device (74), which is a standard 
component of cardiologic echo-doppler systems like the Sonos 2500, Hewlett- 
Packard, USA. 

The information that is transferred to the personal computer (76) for 

10 analysis is: arterial blood pressure measurements (80) from the blood 

pressure monitor (70) including aortic systolic (AoP SYS ) and diastolic (AoP DlA ) 
blood pressures; Volume measurements (82) from the echo-doppler system 
(72) are ventricular end-diastolic (EDV) and end-systolic (ESV) volumes; 
Time measurements are the time T (84) and the time "t^" (86) which are 

15 measured with the echo-doppler system (72). The time T is measured 

between the R-wave of the electrocardiographic signal and the onset of blood 
ejection from the left ventricle. The time "t^" is measured between the R- 
wave of the electrocardiographic signal and the termination of blood ejection 
from the left ventricle. The onset of ejection and the termination of ejection 

20 are determined by the opening and the closing of the aortic valve, 

respectively, which are demonstrated by M-mode echocardiography, or by the 
initiation and termination of a proximal aortic doppler flow waveform. These 
modes are standard on most cardiologic echo-doppler systems. 

These measurements are transferred to the personal computer (76) by 

25 a data transfer module (78) which can be the keyboard of the personal 
computer (76) for manual entry of the parameters, or a standard serial 
communication or network connection for automated data transfer between 
the standard measurement devices (70, 72 and 74) and the computer (76). 
The estimation of the slope and the intercept of the ESPVR from the 

30 non-invasive pressures, volumes and time intervals follows a similar 

procedure as described earlier for the invasive data. The main difference is 
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that now the selection of the time T for the estimation process is limited to 
the region of the onset of ejection. 

Referring now to Figure 7B, the inputs to the personal computer 
analysis program are the systolic and diastolic aortic blood pressures (80), 
the end-diastolic and end-systolic ventricular volumes (82), the time T (84) 
and the time 'W (86). 

The value of the normalized time t N (88) is calculated from Equation 2, 
based on the ratio between T (84) and "t^" (86). The value of the mean 
elastance function at normalized time t N (90) is calculated from the stored 
mean elastance function (92) as detailed above for the first preferred 
embodiment. 

Finally, the slope and intercept of the ESPVR (94) are calculated by 
equations 7 and 8. 

The primary application of this technology is to provide a non-invasive 
cardiovascular assessment tool that can be applied in conjunction to existing 
devices linking echocardiography imaging systems with automated arm-cuff 
and tonometric arterial pressure recording systems. Using pressure-volume 
and flow data that can be assessed non-invasively, the algorithm would 
predict the ESPVR, calculate an effective arterial elastance (E a = P es /SV), 
which is a correlated parameter that characterizes the arterial "load", and link 
these two parameters to predict cardiovascular status and the likely response 
to interventions based on known modeling equations. For example, the slope 
E max , intercept V 0 of the ESPVR, end-diastolic chamber volume EDV, and 
value of arterial load E a , can be combined to predict the stroke volume of the 
heart (Eq. 9), a common parameter that is important in managing patients 
with heart failure. 

SV = E^ x (EDV^AE^+EJ Eq. 9 

Analogous equations have been derived for cardiac work, efficiency, systolic 
pressure, ejection fraction, and other common hemodynamic variables. Each 
of these equations provides the dependent cardiovascular parameter as a 
function of circulating blood volume (EDV), effective arterial load (E a ) which 
also includes heart rate effects, contractility (E max ) and chamber remodeling 
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(V 0 ). Curves describing each of the dependent parameters as a function of a 
varying independent variable can be displayed separately. For example, the 
equation above can be plotted as SV versus EDV, SV versus E ai SV versus 
E^, etc. and for each plot, the exact position of the patient's current status 
5 can be depicted. Figures 8A and 8B illustrates two such curves, where 
variables for cardiac stroke volume SV, and stroke work SW are plotted 
versus arterial loading (E a ) in Figures 8A and 8B, respectively. The current 
status of the patient is noted by the solid circle. The consequence of a 25% 
decline in arterial load, which could be obtained by vasodilator therapy, can 

10 be predicted for both parameters based on this framework. Thus, the 
physician can be shown a series of function curves, which dissect the 
integrated information embodied by the single-beat elastance calculations into 
more practical displays of the cardiovascular system. Repeat curves 
obtained on different office visits could be stored and displayed as temporal 

1 5 histograms to indicate chronic progress of clinical status. This will likely prove 
of great clinical value of the chronic management of heart failure patients. 
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What is claimed is: 

1 . A method of determining an end-systolic pressure-volume 
relationship (ESPVR) of a human heart, comprising steps of: 

measuring a blood pressure of the heart's left ventricle and 
producing a pressure signal indicative thereof; 

measuring a volume of the heart's left ventrical and producing a 
volume signal indicative thereof; 

monitoring the heart to produce ECG timing signals indicative of 

its activity; 

analyzing the pressure, volume and timing signals based on a 
previously determined model of normalized elastance function for a 
normalized time t^, to determine the ESPVR; and 

displaying the ESPVR; 

2. A method according to claim 1 wherein the step of analyzing is 
carried out by a computer. 

3. A method according to claim 2 wherein the pressure, volume and 
timing signals are digital signals. 

4. A method according to claim 2 wherein the pressure, volume and 
timing signals are produced as analog signals and converted to digital signals 
for inputting to the computer. 

5. A method according to claim 1 , wherein the step of analyzing 
comprises the steps of: 

recording the pressure signal and the volume signal over at 
least one cardiac cycle; 

identifying an R-wave of the ECG timing signal; 
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sequencing the pressure signals into pressure signals of 
separate ones of the at least one cardiac cycle based on an identification of 
the R-wave; 

sequencing the volume signals into volume signals of the 
separate ones of the at least one cardiac cycle based on the identification of 
the R-wave; 

group averaging the pressure signals of the separate ones of 
the at least one cardiac cycle to produce a mean pressure; 

group averaging the volume signals of the separate ones of the 
at teast one cardiac cycle to produce a mean volume; 

identifying an occurrence of end-ejection based on one of the 
volume signal from the step of measuring the volume and the pressure signal 
from the step of measuring the blood pressure; 

determining a time interval t max between an occurrence of the R- 
wave and the occurrence of end-ejection; 

calculating an absolute time t based on the normalized time t N 
according to an equation t = t N x t max ; 

determining a pressure at time t r P(t), and a pressure at time 
W . Ptfmax ). based on the mean pressure; 

determining a volume at time t, V(t), and a volume at time t max , 
V(t max ), based on the mean volume; 

determining a normalized pressure P N (t N ) based on the P(t) and 

the Pft^); 

calculating a normalized elastance function for the time t N , 
E N (t N ), based on a stored mean elastance function; 

calculating a volume-axis intercept V 0(SB) of the ESPVR based 
on an equation 

V 0(S B, = [ PN(t N )xV(t max )-V(tN)xE N (t N )] ; and 
[P N (t N ) - E N (t N )] 
calculating a slope E max{SB) according to an equation 

E max (SB) = P(tmax)/[V(tmax) " Vq(SB)]- 
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6. A method according to claim 5, wherein: 

the step of measuring the blood pressure includes measuring 
aortic blood pressure, and 

the step of identifying an occurrence of end ejection is based on 
the blood pressure of the heart's left ventricle and the aortic blood pressure. 



7. A method according to claim 5, wherein: 

the mean elastance function is stored in tabulated form, and 
linear interpolation is used to find E N (t N ) from a nearest 
tabulated point. 

8. A method according to claim 5, wherein: 

a user-controlled value of the normalized time t N , from which 
the absolute time t is calculated, is selected from a range of about 0.25 to 
about 0.35. 



9. A method according to claim 1 , further comprising steps of: 

determining a time interval t between an R-wave of the ECG 

signal and an onset of blood ejection from the heart's left ventricle; 

determining a time interval t max measured between the R-wave 

of the ECG signal and a termination of the blood ejection from the heart's left 

ventricle, 

wherein the step of analyzing the pressure, the volume, and the 
ECG timing signals comprises: 

recording the blood pressure signal, the blood pressure 
signal including aortic systolic A 0 P SYS and diastolic A 0 P D , A blood 
pressures; 

recording the volume signal, the volume signal including 
ventricular end-diastolic EDV and end-systolic ESV volumes; 

calculating a normalized time t N from an equation 

tN ~ t/tnww I 
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calculating a normalized elastance function for the time 
t N , E N (t N ), based on a stored mean elastance function; 

estimating a volume-axis intercept V 0(SB ) of the ESPVR 
based on an equation 

V 0 (SB)-AOP = 

[A 0 P DIA xESV/(0.9xA o P SY s)-EDVx E N (t N )] ; 
[AoPq^O^xAoPsys)- EM\ 
estimating a slope E^sb, based on an equation 

"-max(SB)-AOP 

= (0.9xA o P SYS )/(ESV-V 0(SB) ) ; 
wherein the onset and the termination of the blood ejection are 
determined by one of M-mode echocardiography and by initiation and 
termination of a proximal aortic Doppler flow waveform. 



10. A method according to claim 9, wherein the time interval t is limited 
to a region of the onset of the blood ejection. 

11. A method according to claim 1 , wherein the step of measuring the 
blood pressure is performed by an invasive pressure measurement system. 

12. A method according to claim 1 , wherein the step of measuring the 
volume is performed by an invasive volume measurement system. 

13. A method according to claim 1 , wherein the step of monitoring the 
human heart to produce the ECG signal is performed by an 
electrocardiographic monitoring system. 

14. A method according to claim 1 , wherein the step of measuring the 
blood pressure is performed by a non-invasive blood pressure monitor. 

1 5. A method according to claim 1 , wherein the step of measuring the 
volume is performed by an echo-Doppler system. 



WO 98/19594 



20 



PCT/US97/19695 



16. An apparatus for determining an end-systoiic pressure-volume 
relationship (ESPVR), comprising: 

a blood pressure monitor for measuring a blood pressure of a 
left ventricle of a human heart and producing a pressure signal indicative 
thereof; 

a volume measurement system for measuring a volume of the 
left ventricular chamber and producing a volume signal indicative thereof; 

an electrocardiograph monitor for monitoring the human heart 
and producing an ECG signal of the human heart; 

a computer for receiving the pressure signal, the volume signal, 
and the ECG signal, the computer including: 

means for determining the ESPVR based upon the 
pressure, the volume, and a timing derived from the ECG signal and 
determining the ESPVR. 

17. An apparatus according to claim 16, wherein the determining 
means comprises: 

means for recording the pressure signal and the volume signal 
over a plurality of cardiac cycles; 

means for identifying an R-wave of the ECG timing signal; 

means for sequencing the pressure signals into pressure signals 
of separate ones of the plurality of cardiac cycles based on an identification of 
the R-wave; 

means for sequencing the volume signals into volume signals of 
the separate ones of the plurality of cardiac cycles based on the identification 
of the R-wave; 

means for group averaging the pressure signals of the separate 
ones of the plurality of cardiac cycles to produce a mean pressure; 

means for group averaging the volume signals of the separate 
ones of the plurality of cardiac cycles to produce a mean volume; 
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means for identifying an occurrence of end-ejection based on 
one of the volume signal from the step of measuring the volume and the 
pressure signal from the step of measuring the blood pressure; 

means for determining a time interval t max between an 
occurrence of the R-wave and the occurrence of end-ejection; 

means for calculating an absolute time t based on a normalized 
time t N according to an equation t = t N x t^ ; 

means for determining a pressure at time t, P(t), and a pressure 
at time t^ , ), based on the mean pressure; 

means for determining a volume at time t, V(t), and a volume at 
time t^ , ) t based on the mean volume; 

means for determining a normalized pressure P N (t N ) based on 
the P(t) and the P(t max ); 

means for calculating a normalized elastance function for the 
time t N , E N (t N ), based on a stored mean elastance function; 

means for calculating a volume-axis intercept V 0(SB) of the 
ESPVR based on an equation 

V 0( sb) = Ph(W^(U ' V(t N )xE N (t N )j ; and 
[P N (U - E N (t N )] 

means for calculating a slope E max(SB) according to an equation 

max(SB) 

= P(U)/(V(t max )-v 0(SB) ]. 

18. An apparatus according to claim 16, wherein the determining 
means comprises: 

means for recording the blood pressure signal, the blood 
pressure signal including aortic systolic A 0 P SYS and diastolic A 0 P D1A blood 
pressures; 

means for recording the volume signal, the volume signal 
including ventricular end-diastolic EDV and end-systolic ESV volumes; 

means for calculating a normalized time t N from an equation 
t N = t/tmax • where a time interval t is determined by the volume measurement 
system and the electrocardiograph monitor to be between an R-wave of the 
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ECG signal and an onset of blood ejection from the heart's left ventricle, and 
is input to the computer, a time interval t^ is determined by the volume 
measurement system and the electrocardiograph monitor to be between the 
R-wave of the ECG signal and a termination of the blood ejection from the 
heart's left ventricle, and is input to the computer; 

means for calculating a normalized elastance function for the time t N , 
E N (t N ), based on a stored mean elastance function; 

means for estimating a volume-axis intercept V 0(SB) of the ESPVR 
based on an equation 

Vo( S b).aop = [A 0 P OIA xESV/(0.9xA 0 P 8rs )-EDVx E N (t N )] - 
(A o P D1A /(0.9xA o P SYS ). E N (t N )] ;• 

means for estimating a slope E max(SB) based on an equation 
(O.gxAoPsYs^ESV-V^s^). 

19. An apparatus according to claim 16, wherein: 

the blood pressure monitor comprises an invasive pressure 
measurement system. 

20. An apparatus according to claim 16, wherein: 

the blood pressure monitor comprises a non-invasive pressure 
measurement system. 



21. An apparatus according to claim 16, wherein: 

the volume measurement system comprises an invasive volume 
measurement system. 

22. An apparatus according to claim 16, wherein: 

the volume measurement system comprises a non-invasive 
volume measurement system. 

23. An apparatus according to claim 16, wherein: 
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the volume measurement system comprises a non-invasive echo-Doppler 
system. 
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